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Abstract - We investigate theoretically the dynamics and the charge transport properties of 
a rod-shaped nano-scale rotor, which is driven by a similar mechanism as the nanomechanical 
single- electron transistor (NEMSET). We show that a static electric potential gradient can lead to 
self-excitation of oscillatory or continuous rotational motion. We identify the relevant parameters 
of the device and study the dependence of the dynamics on these parameters. We discuss how 
the dynamics is related to the measured current through the device. Notably, in the oscillatory 
regime, we find a negative differential conductance. The current-voltage characteristics can be 
used to infer details of the surrounding environment which is responsible for damping. 



Introduction. — In recent years it has emerged that 
the coupUng of electrical and mechanical degrees of free- 
dom on the nanometer scale provides the opportunity to 
build novel devices, extending the concepts of conventional 
electronic devices [Ull] . A seminal example is the nanome- 
chanical sing Ic-clectron transistor (NEMSET) [SHS], where 
electrons can tunnel from a source to a drain electrode via 
a movable island or grain, whose charge (i.e. the number 
of electrons that occupy the island) is determined by the 
Coulomb-blockade effect. Since the tunneling amplitude 
depends exponentially on the position of the island, the 
current is very sensitive to the mechanical motion. For a 
sufficiently large bias voltage a self-excitation of periodic 
oscillations occurs in conjunction with charging and de- 
charging of the island. This mechanically assisted charge 
transport is called electron shuttling [3] . The grain is em- 
bedded in a medium which creates a restoring-force and 
determines essentially the eigenfrequency of the oscillation 
of the shuttle. 

In this letter we consider a device which can rotate freely 
and which is driven by the same mechanism as the elec- 
tron shuttle described above. The driving force is deter- 
mined by a static voltage, which is applied along the device 
[THS]. In contrast to the conventional shuttles this results 
in a force that depends non-lincarly on the relevant sys- 
tem coordinate. The coupling of mechanical motion and 
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tunneling leads, as we will show, to the self-excitation of 
oscillatory and rotational motion even in the presence of 
damping. The frequency of oscillation (rotation) depends 
on the ratio of the driving force and the friction. 

While being based on the same principle as the charge 
shuttle, the present device exhibits markable differences. 
In the oscillatory regime, the current through the device 
decreases with increasing bias voltage, which is surpris- 
ing from the conventional shuttling point of view. We 
attribute this effect to the presence of a separatrix in the 
phase space, which separates oscillations from rotations. 
Approaching the separatrix involves a slowing down of the 
dynamics and decreasing oscillation amplitudes, which re- 
sults in the decreasing current. In the rotational regime, 
one may realize a nano-scale motor, which is driven by a 
static voltage [9]. 

Using numerical calculations supplemented by an an- 
alytical analysis, we find that the dynamics of the rotor 
is governed by three dimensionless parameters: tunneling 
length, field strength and damping constant. This circum- 
stance allows us to predict, for example, a transition from 
oscillatory to rotational motion fairly independent of the 
actual realization of the device. Our analysis provides an 
intuitive picture of the dynamics for a large range of pa- 
rameters. 

Model. — A possible realization of the device, dis- 
cussed in the present letter, is sketched in fig.[T] It consists 
of two charge islands (e.g. quantum dots, metallic nano- 



p-1 



A. Croy and A. Eisfeld 




Fig. 1: (Color on-line) a) Artists view of a possible realiza- 
tion of the single electron motor using (multi-wall) carbon nan- 
otubes and two charge islands, b) Sketch of the relevant parts 
and coordinates. 

clusters, etc.), which wc denote by A and B, placed at the 
ends of a rigid rod of length i. This structure can rotate 
in the x-y plane around an axis in the middle of the rod. 
The support shaft of the rotor could be realized similar 
to the one of ref. |10| by multi-walled carbon nanotubes 
or in a way proposed in ref. [7]. The motion of the rotor 
around the pivot axis is characterized by the angle Q and 
the corresponding moment of inertia /. 

To drive the device, the rotor is placed between source 
(S) and drain (D) contacts which are connected to elec- 
tron reservoirs. The contacts arc kept at different chemi- 
cal potentials by an externally applied bias voltage. The 
electrons may tunnel between the contacts and the charge 
islands. The respective tunneling amplitudes depend sen- 
sitively on the position of the islands. Specifically, we 
assume an exponential dependence [3]. For example, 
the tunneling amplitude from the source to island A is 
Ts^A exp[— Axs.a/As.a], where Axs.a is the distance 
between the source contact and island A and As. a is the 
tunneling length. Similar expressions are used for the tun- 
neling amplitudes Td.a, Ts.b and Td,b- It is convenient 
to introduce the dimensionlcss ratio 

^S.A = , (1) 

which can be regarded as a measure of the change of dis- 
tance with respect to the tunneling length. 

Since the islands are nanoscale objects, their charging 
will be determined by the Coulomb-blockade effect [TTJ 



For sufficiently small tunnel coupling electrons are 
transferred sequentially. According to the orthodox theory 
of Coulomb blockade [TT]|13II14| the tunneling rates are (for 
island A) 

W^s,A(e)=rsAe-«^'*™^(®) , (2a) 

w^D,A(e) =rDAe+«° *™^(®) , (2b) 

where for example Fsa is defined as the tunneling rate 
from source to dot A when the rotor is in the perpendicular 
position, i.e., for 8 = 7r/2. Similar relations hold for the 
tunneling rates of island B (replacing cos(0) by cos(0 -|- 

In the following, we consider a very large charging en- 
ergy such that only one excess electron can occupy each 
island at the same time (Coulomb blockade regime). Fur- 
ther, for simplicity, we assume that the bias voltage is 
very large (large bias limit). The large bias approximation 
holds as long as the respective broadened energy level is 
completely in the transport window, i.e., for source-drain 
voltages eVsD > k^T and eVsD > Fsa/b + Tda/b- Here, 
A:b is the Boltzmann constant and T denotes the tempera- 
ture of the electronic reservoirs. Then, electron transport 
occurs only in one direction [151 : From the source contact 
electrons can tunnel onto the islands. On the other side 
they can tunnel from the islands to the drain contact. 

Dynamics. — Denoting the electronic population of 
island A (B) by -Pa(b) wc can write for average charges 
of the islands (7a(B) = ^ePA(B) where Pa(^) can attain 
values between and 1. Due to the tunnel coupling to 
the leads the populations Pa{B) time-dependent and 
determined by the rate equation 

PA = Ws,Am{l-PA)-WD,A{^)PA- (3) 

The first term on the right hand side describes tunnel- 
ing from the source onto island A and the second term 
is responsible for tunneling from the island to the drain 
contact. 

When island A possesses a charge qa then it experiences 
a force in x-direction which is proportional to the elec- 
tric field strength £ ex Vsd induced by the source-drain 
voltage Vsd between the contacts: — qa^- Then, in 
the mean field description introduced above, the torque 
acting on the rotor is given by Ma = —e Pa £ § sin(9). A 
similar expression holds for island B (replacing sin(0) by 
sin(e -I- tt)). 

To keep the discussion transparent, in the following we 
will take all tunneling rates to be equaQ and denote them 
by r. This rate can then be used to introduce a dimen- 
sionless time, 

T = Tt, (4) 

^ In practice the couplings of the different islands to the two 
leads will not be perfectly identical. From our numerical simulations 
we have found that the features discussed in the present work are 
quite robust with respect to changes in the relative couplings. For 
example, even for very asymmetric tunneling rates T a,S A,D ~ 4 
the results presented in the following, remain qualitatively the same. 
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Fig. 2: (Color online) Dynamics of the rotor for various rj and fixed 7 — 0.01 and ^ = 2. (a) Phase portrait of the steady-state 
dynamics. Within the gray-shaded area an oscillatory motion takes place. The values of 77 are given next to the respective 
trajectory, (b) Exemplary time evolution that leads to the corresponding steady-state dynamics, (c) Sketch of the dynamics 
for the four cases shown in b). 



and a dimensionless field strength, 
77 = e£i/{2ir^) 



(5) 



Due to the torque Ma + Mb there will be a change in 
angular momentum L of the rotor. For the dimensionless 
angular momentum 

A - — (6) 
IT ^ ' 

the equation of motion reads 

= -77 (sin(e)PA + sin(e + ^)Pb) - F{K) 
dr (7) 

= -?7sin(e)AP-F(A) , 

where we have added a phenomenological damping term 
F{A) and in the last line we have defined AP = Pa — Pb- 
In the following calculations we choose the damping to be 
linearly dependent on the velocity [71ll6j. i.e., F{A) = 7A, 
with a dimensionless damping constant 7. The damping 
will typically be accompanied by fluctuations, whose in- 
fluence will be discussed in a later section. 

The population difference AP, which appears in eq. (O, 
also depends on time. We have 



dAP 



where we have also set all tunneling lengths equal to ^. 
Finally, the angle Q is related to the momentum A by 



a7 



A 



(9) 



-2 sinh[^ cos(e)] - 2 cosh[^ cos(e)] AP , (8) 



The three coupled equations, ([7]), © and govern the 
dynamics of the system. In the following, we will investi- 
gate how this dynamics depends on the parameters ^. 7 
and rj. We are mainly interested in the steady-state dy- 
namics at long times after the initial irregularities have 
vanished. 

For given damping 7 and tunneling length ^ we find 
three qualitative different regimes of the steady-state mo- 
tion depending on the values of the driving 77. These are 
illustrated in fig.Ht. In fig.[5^,b numerical solutions of the 
set of equations ©-([S]) are shown for 7 = 0.01 and C = 2. 

In fig. the steady-state solutions for various values of 
77 are depicted by representing the angular momentum A 
as a function of the angle Q. Different initial conditions 
lead to the same steady-state solution up to the direction 
of rotation, which is opposite for and — 80- In fig. 
^Bp examples of the initial time dependence of the angular 
momentum is shown for a situation where the rotor is 
initially at rest at an angle 9 = 0.4. 
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Fig. 3: (Color online) Period of steady-state oscilla- 
tions/rotations is shown as a function of 77/7. Different symbols 
(colors) indicate different values of 7, as given in the legend. 
The inset shows the corresponding maximal angle relative to 
e = 7r/2. 

For 77 < 7 the rotor shows a relaxation towards the 
equilibrium position given hy Q — 7r/2, i.e., both islands 
A and B are in the middle between the source and the 
drain contact. Obviously, in this regime the dynamics of 
the rotor is dominated by the damping, which leads to the 
motion of the rotor being slow compared to the tunneling 
time. Consequently, the populations are in equilibrium 
with the reservoirs at each instant of time. Additionally, 
the strong damping leads to decreasing amplitudes and, 
in the symmetric case considered here, the population dif- 
ference vanishes accordingly. Since the rotor is not ac- 
celerated anymore, the motion eventually grinds to a halt. 
Electrons will tunnel with the constant rate F to/from the 
islands to the contacts. The stationary solution is given 
by the single (blue) point in fig. [2^, located at Q = tt/2 
and A = 0. 

Upon increasing i] one reaches a regime where the rotor 
oscillates around the position 6 = 7r/2 and the force on 
the charged islands is not yet strong enough to lead to a 
rotational motion. This regime is characterized by closed 
curves within the shaded area in fig. ^j^. The outmost tra- 
jectory of this area was obtained from an adiabatic approx- 
imation to eq. ([5]) and considering extremal swing. The 
occurrence of the oscillatory regime is a direct consequence 
of the distance dependence in the tunneling rates. For suf- 
ficiently weak damping, inertia leads to slightly increasing 
amplitudes in each period. Due to the increasing tunnel- 
ing rates, this comes along with an increase of the popu- 
lation difference and therefore a larger backward acceler- 
ation. Of course the increase of the oscillation amplitude 
will be limited by the damping, which becomes more and 
more important as the velocity (or amplitude) increases. 
Consequently, one has two competing mechanisms: one 
that favors to increase the oscillation amplitude (electric 
field/tunneling) and one that acts to decrease it (damp- 
ing). The result is the occurrence of stable oscillations 
in the long-time limit. Notice, that these oscillations are 



akin to the self-excited shuttling in NEMSET [3]. How- 
ever, in contrast to the conventional charge shuttle, for the 
rotor the period of the oscillation increases with increasing 
77. This behavior is illustrated in fig. [31 where the period 
of steady-state oscillation/rotation is shown as a function 
of 77/7. In the oscillatory regime, indicated by the gray- 
shaded clear increase of the period of oscillation 
can be seen, spiking at 77/7 « 2.8. This point corresponds 
to the separatrix as shown in fig. [5^. 

For even larger values of 77 periodic rotational motion 
sets in (see the curves for = 0.04 and = 0.05 in fig. [2^). 
The direction of rotation in the long-time limit depends 
on the initial conditions. In this regime, the motion is 
dominated by the driving through the electric field. The 
rotational energy becomes large enough to overcome the 
damping and to allow an almost free rotation, character- 
ized by a constant angular momentum. This case is shown 
in fig. [2^ in the upper panel. 

Angular momentum. — With respect to a possible 
application as a motor, the temporally averaged angular 
momentum in the steady-state is an important quantity. 
Thus, in our numerical simulations we calculate 

(A) = - A(t)dt , (10) 

where we have chosen the time ti large enough to ensure 
that we are in the regime of stationary dynamics and we 
have taken the time interval T to cover many periods of 
oscillation (rotation). In fig. |4] numerical results for var- 
ious ^ and 7 are shown. The different colors of the dots 
represent different values of 7 for the same ^ as in fig. O 
As can be anticipated from eq. ([7]), with F{A) — 7A, the 
curves for different values of 7 nearly coincide. 

In the upper row of fig.|l]the dependence of |(A)| on 7^/7 
is shown. Obviously, in the damping dominated regime 
the average angular momentum is zero, since A in 
the stationary limit. Also in the oscillatory regime (in- 
dicated by the gray-shaded area) one has (A) = 0, since 
the angular momentum changes its sign in each period 
(see, e.g., the third row of fig. [21d). Thus, the averaged 
A is not an ideal quantity to gain insight into the oscilla- 
tory regime. However, it is well suited in the rotational 
regime, which is characterized by a finite value of the av- 
erage angular momentum (A) and its modulufH. In the 
left column of fig. 2] the focus is on relatively small values 
77/7 < 10, in the right column the region from 77/7 = 10 
to 77/7 = 150 is covered. In this region the dependence 
of the angular momentum on 77/7 can be estimated an- 
alytically: In the stationary limit, the time-averages of 
the dissipated work and the work due to accelerating the 
charges have to be equal, (Wfrk) = (Wchargc)- Chang- 
ing the time-average to an angular average amounts to 

^Whcn considering (A) instead of |(A)| a branching of the curve 
appears, where the upper and lower branches describe anti-clockwise 
or clockwise rotation, respectively. The direction of rotation depends 
in a non-trivial way on the initial angle ©o- Due to the symmetry 
of system, changing 0o to — 0o reverts the direction. 
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7 /p de A(e) = -77 /q dQ AP(e) sin (9). For large val- 
ues of 77/7, the angular momentum is constant, (A) w 
A = Agtat- After integration by parts and replacing the 
resulting derivative of the population difference by eq. ([S]) , 
dAP/de = {dP/dT)/{dQ/dT), one obtains 



A.tat«,/c(C);^-d(0 



(11) 



with c(0 = (cos (9) 2sinh(^cos(9)))e and d{^) = 
(cos (9) 2cosh(^cos(9)) AP)q. Using the Jacobi- Anger 
identity, the first term becomes c(^) = 27i(^), where /„ is 
the nth order modified Bessel function of the first kind. 
The second term is more difficult to evaluate. To elimi- 
nate the population difference, one has to repeatedly use 
eq. (|5]). However, in the limit 77/7 — >■ cx) this term vanishes. 

We have found that these analytical estimates describe 
our numerical results very well in the region of large 77/7. 
For the case ^ = 2, shown in figs. [2] and |4] we have numer- 
ically found c(0 ~ 3.18 and d{^) w 19.76. We found also 
very good agreement for other values of ^. 

Current. — Since the dynamics of the rotor is essen- 
tially driven by charge transport it is instructive to take 
a closer look at the electric current. The time-dependent 
current (in units of the electron charge e) is given by 

Jit)/e = M/A,D(e(t)) PA{t) + WB,u{Oit)) Psit) . (12) 

As before we are primarily interested in the time-averaged 
current (J), which we define in the same manner as in 
eq. (fTOj) . Note, that the resulting stationary current is 
independent of the initial condition and, in particular, of 
the direction of rotation. 

Numerical results for ^ = 2 are shown in the second row 
of fig. m Here the average current is plotted as a function 
of 77/7 for different 7 (the same as for the angular momen- 
tum in the upper row). Again, all curves coincide quite 
well. From the behavior of the averaged current one can 
clearly distinguish the regime where the rotor oscillates 
back and forth (corresponding to the shaded area) and 
the regime corresponding to a motionless rotor (7^/7 < 1). 

For a very small ratio 7^/7 the current is constant, since 
the rotor does not move and electron transport takes place 
by tunneling to and from the islands in the middle between 
the contacts. In contrast, in the region where the rotor os- 
cillates one finds a rapid decrease of the average current. 
When the rotational motion starts the current rises again 
slowly. The decrease in the current is related to the fact, 
that the mechanical motion is becoming slower as the sep- 
aratrix, which separates the oscillation trajectories from 
the rotating ones, is approached (cf. fig. With increas- 
ing T//7 the amplitude of the osciUations increases and 
the rotor spends an increasing time at the reversal point, 
where the tunneling to the opposite electrode is strongly 
diminished. This results in a decrease of the current. 

In the limit of a very large ratio 77/7, the stationary 
current attains a constant value, which is readily calcu- 
lated from eqs. (O and (|12p by replacing the rates by their 
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Fig. 4: (Color online) Dependence of the average angular mo- 
mentum (A) and the average current (J) on the ratio 77/7. 
Symbols denote results for ^ = 2 and different values for 7 as 
given in the legend. Right column: focus on large values of 77/7. 
Left column: details for small 77/7. The dashed line indicates 
the behavior obtained from eq. and the dotted line shows 
the asymptotic value of the stationary current ( J)stat = F /o(C)- 
The inset shows (J) for the same 77/7, but on a larger scale. 



time averages. One obtains (M^s,a) = Fsa7o(Cs,a) and 
(W^d,a) = rDA7o(^D,A), where /o(0 denotes the modified 
Bessel function of the first kind. The stationary current 
(for the case of identical tunnel parameters) is then given 
by (J)stat = r^o(0- This value is depicted in fig.|4]as the 
horizontal dotted line. It is in very good agreement with 
the numerical results. 

In addition to the case ^ = 2, in fig. 2] also numerical 
results for other values of ^ are shown. One sees that for 
larger ^ the drop in the current sets in for smaller values 
of 77/7 and is more pronounced. 

Robustness of the results. — In the following, we 
would like to briefly comment on the robustness of our 
results with respect to asymmetries in the couplings to 
the leads and fluctuations of the environment (e.g. related 
to internal friction or interaction with some background 
gas at finite temperature). In the present work we have 
restricted ourselves to an ideal situation, neglecting these 
effects, since we wanted to keep our analysis as clear as 
possible. Our results are only weakly affected by includ- 
ing different couplings and fluctuations. For the couplings 
this has been discussed in connection with eq. To 
gain insight into the influence of fluctuations we included 
a Langevin term in eq. ([7]) and solved the resulting stochas- 
tic differential equation numerically. For small 77 the fluc- 
tuations have the largest effect. In the oscillatory regime 
and close to the separatrix they can lead to sporadic rota- 
tions. Similarly, in the standstill regime oscillations with a 
few periods can be triggered, resulting in a small decrease 
of the current. Overall, this leads to a smoothing of the 
observed features for the time-averaged quantities. How- 
ever, for moderate fluctuation strengths (approximately 
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smaller than the dissipative energy 7(A)) the qualitative 
behavior is not changed. In the rotational regime (large 77) 
our results were (as expected) nearly unaffected by fluc- 
tuations. 

Conclusions. — We have investigated the dynamics of 
a nanoscale rotor which is driven by the sequential tunnel- 
ing of electrons between electronic contacts and the device. 
The interplay of the electric and the mechanical degrees of 
freedom leads to three different regimes for the motion of 
the rotor, which can be obtained by adjusting the external 
bias voltage: stationary, oscillatory and rotation. 

The oscillatory regime has some similarities with the os- 
cillatory motion of the conventional charge shuttles. How- 
ever, in contrast to the charge shuttles, the motion of the 
rotator is slowed down with increasing driving strength. 
This behavior results in a decrease of the electric current 
with increasing bias voltage. In the rotational regime the 
current starts to increase again and the device can be con- 
sidered as a nano-scale motor. It is worth noting that the 
setup discussed in the present paper has some similarities 
with bio- motors [8], such as the FO motor of adenosine 
triphosphate synthase. Further examples of molecular ro- 
tors can be found in ref. |17j . 

One might speculate about possible applications of the 
presented device. Due to the high sensitivity of the current 
to the ratio 77/7, the device may be used to detect changes 
in its environment. For example, the damping constant 7 
may be determined by tuning the voltage and maintaining 
the same current. In this scenario the rotor can be used 
as a sensor. Another possible application is given in the 
context of electron pumping. Especially in the rotation 
regime the charge transport is determined by the mechan- 
ical motion. This provides the possibility of a mechan- 
ically stabilized pumping in the high frequency regime. 
In view of recent achievements regarding the cooling of 
nanomechanical systems to the quantum regime [THHS], 
the questions if and how it is possible to cool a nanome- 
chanical rotor and its interplay with quantum gears |22| 
will be an interesting subject for further studies. 

Altogether, the investigated nano-scale rotor is an in- 
teresting nanoelcctromechanical device with potential ap- 
plications as a sensor or a single-electron motor. 

* * * 
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